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Abstract 

For a one- dimensional model in which the two-body interactions are long- 
range and strong, the system almost crystallizes. The harmonic modes of 
such a lattice can be used to compute the ground state wave function and 
the dynamical density-density correlations. We use this method to calculate 
the density corrrelations in the Calogero-Sutherland model. We show numer- 
ically that the correlations obtained are quite accurate even if the coupling 
is not very large. Our method is considerably simpler than the ones used to 
derive the exact results, and yields expressions for the correlations which are 
easily plotted. Our equal-time correlations can be expanded in powers of the 
inverse coupling; we show that the two leading order terms agree with the 
exact results which are known for integer values of the coupling. 
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I. INTRODUCTION 



The Calogero-Sutherland-Moser model (CSM) [1-3] has received much 
attention recently due to its connection with a variety of interesting problems. 
Some examples are random matrix theory [Q, quantum spin chains with 
long-range interactions Luttinger liquids 0], Gaussian conformal field 

theories M, edge states in a quantum Hall system H, generalized exclusion 



statistics [10-16], and nonlinear internal waves in a stratified fiuid [0. 

Although the quantum mechanical energy spectrum of the CSM has been 
known for a long time [1-3], its collective properties [18-22] and dynamical 
density-density correlations [23-28] are not yet expressed in an easily acces- 
sible form. The correlations are best understood for three particular values 
of the coupling parameter A = 1/2, 1 and 2 which are closely related to 
random matrix theory Q . An impressive amount of information is also 
available for all integer [25-27] and rational [|5H| values of A. However, even 
though these are all exact results, the mathematical expressions are not very 
convenient for plotting, say, the space-time dependence of the correlations, 
unless A (or its inverse) is equal to a fairly small integer. 

It therefore seems useful to pursue simple and approximate methods for 
obtaining the correlations for general values of A. In Sec. II, we present a 
method for doing this which is actually exact in the limit of large A. In that 
limit, the particles almost freeze into a lattice [P, |18]; the harmonic modes 
about such a lattice may then be used to compute the ground state energy, 
wave function and correlations. We use this idea to obtain the first two 
terms in an expansion in 1/A for the ground state energy and the asymptotic 
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expression for the equal-time density-density correlation; these agree with 
the known exact results . We show that the equal-time correlations given 
by our method agree reasonably well with the exact result for all distances 
even if A is not very large. The djTiamical correlations are also obtained with 
almost equal ease by our method, and are explicitly given for any space-time 
point. We discuss the low-temperature specific heat, and make a suggestion 
for improving the correlations by using a slightly modified phonon dispersion. 
In Sec. Ill, we briefly consider a more general model in which the two-body 
interaction decays as the power a of the distance. We find that the system 
behaves as a Luttinger liquid if a > 1, and as a Wigner crystal if a < 1; the 
Coulomb problem with a = 1 lies on the border between the two possibilities 



II. CALOGERO-SUTHERLAND-MOSER MODEL 

We will consider the form of the CSM in which particles on an infinite line 
interact pairwise through an inverse-square potential. The model may also 
be defined on a circle with periodic boundary condition 0; the two versions 
of the model have identical physical properties in the thermodynamic limit 
in which the the number of particles and the length L of the line (or circle) 
are simultaneously taken to infinity keeping the particle density po = Njl^ 
fixed. The Hamiltonian for particles on a line is given by 

where the dimensionless coupling A > 0. To make the problem well-defined 
quantum mechanically, we have to add the condition that the wave functions 



go to zero as \xi — Xn\^ whenever two particles / and n approach each other. 
For A = and 1, the model describes free bosons and free fermions respec- 
tively. Since the two-body potential is singular enough to prevent particles 
from crossing each other, we can choose the wave functions to be either sym- 
metric (bosonic) or antisymmetric (fermionic). The energy spectrum is of 
course the same in the two descriptions. Two of the exactly known results 
for this model are as follows 0, 0]. In the thermodynamic limit, the ground 
state energy Eq is given by 

Eo ^ n'n'X^pl 

N 6m ' ^ ^ 

The sound velocity v is given by the group velocity of long- wavelength modes 
(wavelengths much bigger than the average particle spacing 1/po); 



m 



(3) 



A. Ground State Energy and Wave Function 

Let us now consider the limit in which the coupling A — oo. Then the 
potential energy term in (|l|) dominates over the kinetic term. The potential 
energy is minimized if the particles sit at the sites of a lattice, so that the 
mean position of the n-th particle is (x^) = nj |TB[. The total energy 
at this order in A is then given by 

£o ^ 7r^^^A(A - \)pI 

N 6m ■ ^ ^ 

To the next order in A, we have to consider the harmonic oscillations about 
the mean positions. If we denote the fluctuation by i]n = Xn — n/po, the 



Hamiltonian for the phonons is obtained by Taylor expanding the potential 
energy in (|lD to quadratic order. Thus 

_ pI 3h'X{\ - 1)pI (r/z - Vn? ... 

" V 2m + m t i^-^y ' ^ ^ 

where we have ignored higher order anharmonic terms. We now find the 
dispersion relation ujq for phonons with wavenumber q to be 



UJ 



2 12ti^X{X-l)p^ ^ l-cos(gn/po) 



E 



= — —2 — [ ^/^o|g| - ^ j ' kl < ^Po • (6) 

(Note that q can vary from — vrpo to vrpo in units of 271/ L, so that the total 
number of modes is N). The sound velocity is given by f = {dujq/dq)q=Q = 



Ti'^Poy A(A — l)/m. We see that this agrees with Eq. (||) to leading order in 
A; we will therefore simply replace A(A — 1) by A^ in the expression (|^) for uj^. 
Namely, we use the exactly known sound velocity to correct the coefficient 
of our phonon dispersion. Henceforth we take 

= — ( Trpo |g| - |- ) • (7) 
The zero point energy of the phonons is given by 

= L — ^ (8) 

g 2 J—npo 2tt 2 

When this is added to (^), we precisely recover Eq. (|]); thus the harmonic 
lattice approximation gives the ground state energy correctly upto order A^ 
and A. To get the next order term in the energy, i.e. order 1, we would have 
to use the original dispersion as well as consider some of the anharmonic 
terms beyond Hph, we will not pursue this here. 



We now use the harmonic approximation to write down the ground state 
wave function following Ref. [Q. The unnormalized wave function for a 
collection of decoupled phonons is given by 

^0 = exp ( - ^ ^ UJgrjqTJ.q ) , 

oo 

V, = Vn e'"^'"' . (9) 

n=— oo 

Using the dispersion (|^), we find that 

*o = exp ( - — —--^ ) . (10) 



2 t^n (^-^) 

We see that this is equivalent to the exact ground state wave function of the 
Hamiltonian (|1|) 



^0 = n \P0{XI-Xn)\^ , (11) 
l<n 



if we Taylor expand 



\og{po\xi-Xn\) = log \l-n\ + —f—^ 2{l-ny " ' ' ^ ^ 

Before ending this subsection, let us write down the second quantized 
form for the Heisenberg operators x„(t), 



, ^ n f^po dq ^ 



i{ — —LU„t) , t —i( — —LO„t) 

ttg e^po ^ ' + ajj e ^"o " ' 



Here 

[ a, , oj, ] = 27T6{q - q') . (14) 

We can verify that Pn(t) = mdxnit) / dt satisfies the equal-time commutation 
relation = iMin- 
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B. Dynamical Correlation Function 



We now use the harmonic lattice approximation to derive the dynamical 
density-density correlations defined as 



g{x,t) = (0|p(x,t)p(0,0)|0) 

p(x,t) = ^ - Xn{t)) , 



(15) 



where Xn{t) is the Heisenberg operator given in ([l3|) . To compute this, we 
write the 5-functions in (O) as 



5{x - Xnit)) 



(16) 



'-00 2tx 

We then evaluate expectation values of the form (0| exp i{qxn(t) — q'xi{0))\0) 
using the Baker-Campbell-Hausdorff formula, 



^[B,A]/2 



(17) 



if [B, A] commutes with both A and B. Due to the logarithmic divergence 
at small momenta of an integral of the form / dk/uk, we find that this ex- 
pectation value vanishes unless q = q' . We finally obtain 



g{x,t) 
where 



Po 



F{n,t) 



dq V n q^ 

— exp iq{x-—) - -^F{n,t) 



Po 



Po 



'^po dk hpo 



1 _ e— ^* cosf— ) 



(18) 



(19) 



/o IT 2muJk 

It can be shown that the real part of F{n, t) is always positive, unless both 
n and t are zero. We can therefore perform the q integration in ([18|) to get 



g{x,t) = pI 



4nF{n,t) 



1/2 



exp 



(xpo - nf 
AF{n,t) 



(20) 



The expansion in (|T^) can be understood in a physically intuitive way as 
follows. Each lattice point n contributes a Gaussian (correctly normalized 
to unity) to g{x,t); the squared width of the Gaussian is given by the mean 
square fluctuation 

2F{n,t) / . . n ^2 



Pi 



( Xn{t) - Xo{0) 



Po 



|0) 



(21) 



Let us now study the equal-time correlations g{x, 0) = g{x) in some detail. 



From (p!9|), we see that F{n, 0) is purely real, and -F(0, 0) = 0. We therefore 
obtain 



g{x) = Po S{x) + Po E ( 
where 



n,0) 



1/2 



exp 



(xpo - 
4F(n,0) 



(22) 



AF(n,0) 



27r2 Jo 
This gives us pT | 

AF(n,0) 



dy 



1 — cos (yn) 
y - 1/V2vr 



27r2 7o 



, 1 — COS (yn) 
dy 

y 



(23) 



1 



27r2 



log (27re'^n) - Ci (27rn) 



(24) 



where 7 0.57722 is Euler's constant, and Ci denotes the cosine integral. 
For large integer n, (03) has an asymptotic expansion beginning as 



AF(n,0) 



where 



2n 



C 



1 1 

log [Cn) + 



2'Ke' 



11.191 . 



(25) 



(26) 



The expansion in (PSQ converges extremely rapidly; for = 1 and 2, the 
difference between (p¥) and the first two terms in (pS]) is only about 0.1% 



and 0.01% respectively. For numerical purposes, therefore, we will use (p5|) 
to compute g{x) from (p2D . 

We show a plot of g{x) / pi vs. xpo for A = 2 in Fig. 1; the solid line denotes 
the exact result H], while the dotted line shows our lattice approximation 



(1221) . The agreement appears to be reasonable even though A is not very 



large; for instance, the root mean square fluctuation for nearest neighbor 



particles is given by ( pT]) and (p5|) to be y0.247/A/po, which is as large as 
35% of the lattice spacing l/po for A = 2. In Fig. 2, we show the lattice 
approximation for g{x)/pl for A = 4. It is clear that the system shows an 
increasing tendency to crystallize as A increases. We should note that the 
5-function at the origin has not been shown in Figs. 1 and 2; thus the g{x) 
shown in the figures integrates to one hole, 

/■n/2 

jirn 2 j dx g{x) = po {n - 1) . (27) 

We also note that in the lattice approximation, g{x) does not quite vanish at 
a; = as it should in an exact calculation; however g{0) does become small 
rapidly as A increases. 

We can Fourier transform g{x) to obtain the static form factor 

/oo 
dx e-^«^ g{x) . (28) 
-oo 

Since g{x) 1 for large x, S{q) has a 5-function at the origin. Hence it is 
easier to compute 



S{q) - 271 pI (5(g) = po + Po dx cos (gx) [ ^ 



4F(n,0) 



nj^o JnF{n,0) 



- 2 
(29) 



We show a plot of S{q)/po vs. q/po in Fig. 3 for A = 2 (the sohd and 
dotted hnes again denote the exact result and the lattice approximation 
respectively), and in Fig. 4 for A = 4. We will comment on the divergences at 
g/po = 27r in subsection C. (The small wiggles in Fig. 4 are due to numerical 
inaccuracies; the integral in (^) converges very slowly if A is large. The 
power-law which governs the convergence will also be derived in the next 
subsection) . 

We now consider the dynamical correlations g{x,t). It is convenient to 
introduce a dimensionless time variable 



i = vtpo . (30) 



Then 



1 — cos (yn) exp —itiv — ■u^/27r) 
AF(n,t) = — / dy (31) 



271^ Jo y — y'^/2'K 

In the lattice approximation, plotting g{x,t) is almost as easy as plotting 
g{x). In Figs. 5 and 6, we show g{x,t)/p2 vs. xpo for i = 0.2 and i = 2 
respectively, both for A = 2. The solid and dotted lines show the real and 
imaginary parts respectively. Note that the 5-function at the origin which 
was not shown for g{x) in Fig. 1 has now spread out and become visible at 
finite time in Figs. 5 and 6; hence there is no hole in the real part of g{x, t). 
For a small value of time t and for x close to the origin (Fig. 5), the 5- function 
at the origin spreads out as a free particle with dispersion Uq = hq'^/2m; this 
explains the large oscillations in both the real and imaginary parts. For a 
larger value of time (Fig. 6), and also for larger values of x for any time, the 
particle feels the harmonic force of the lattice. Hence the behavior of g{x,t) 
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is different from the spreading of a free particle; in particular, the large 
oscillations get damped out. Finally, note that for t = 2, that g{x,t) can 
differ appreciably from g{x, 0) only in the region xpo = to 2 since phonons 
can only propagate that far in that time. This is most strikingly visible in 
Fig. 5 in the imaginary part of g{x, t) which vanishes rapidly beyond xpo = 2. 

Given g{x,t), we can Fourier transform the x coordinate to obtain a 
function which we denote by 

oo • 2 

Siq,t) = Po E exp [ - ^ - ] . (32) 

n=-oo Po Po 



Thus the function S{q) given in ( pSf ) is equal to S{q, t = 0). If we now Fourier 
transform (R^) in time, we obtain the dynamical structure factor 



S{q,uj) = / dt e*^* S{q,t) . (33) 

J —oo 

This quantity can be written in terms of all the excited states |z/) of the 
system, 

S{q,u;) = Po E 2n6{u;-E, + Eo) \{u\p{q)\0)\\ 

p{q) = f dx p{x) e-^'^" = E e"''"" • (34) 



It is useful to consider the moments of S{q,uo) p3| 

r°° duj 



/oo Hf ) 
— oo-Siq^u). (35) 
-oo /vr 



These can be numerically obtained in a straightforward way from (|32[) , since 
In{q) = i"{d"'S{q, t)/dt"')t=o- We have already considered Io{q) = S{q) above. 
Since F(0,0) = and {dF{n,t)/dt)t=o = iSn,o ^Po/2?^, we find that the first 
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moment 

IM = (36) 

is independent of A. This is an exact result following from the velocity 
independence of the two-body interactions, and it is called the /-sum rule. 
The second moment l2{(l) is shown in Fig. 6 for A = 2; we have actually 
plotted l2{q)m? / {Ti^ plq^) which is dimensionless and goes to a constant at 
g = 0. 

At this point, we observe the phenomenon of saturation by sound modes 
for small momenta P3|. Namely, for q « 2npQ, S{q,uj) is dominated by 



values of u close to the phonon energy Uq given by (^; thus /„(g)//o(g) 
approaches as q 0. Hence the ?/-coordinate in Fig. 6 should approach 
'7rA/2 at g = which it does. For the same reason, the curves in Figs. 3 and 
4 are linear near the origin with a slope of 1/ (27rA), since S{q)po approaches 
7iq'^/{2mujq). 

The saturation by phonons at low energies also has a bearing on low- 
temperature thermodynamic properties of the CSM like the specific heat 
and pressure. We use the dispersion (0) to compute the free energy per unit 
length / taking the phonons to have zero chemical potential. Thus 

/?/ = r ^\n{l - e-^^^ ) , (37) 
j-oo zvr 

where (3 = l/ksT. After evaluating this, we obtain the specific heat per unit 

length Cv = —Td^f/dT^ to second order in T. We find 

Cy = !^ + 6C(3)A fc|T2 
3hv 7T hmv^ ' 

C(3) = E - • (38) 
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On comparing this with the exact result in Ref . |T5[ , we see that we have the 
correct coefficient for the term of order T; this is related to having the right 
phonon velocity at g = 0. However we have a coefficient A instead of A — 1 
for the term of order in (pH]). A different approach to the GSM based on 
collective field theory actually gives a dispersion 



-(^PoA|g| - (A-1)^), (39) 



which leads to the correct coefficient A — 1 for the term. All this suggests 
that we may get more accurate low-momentum or long-distance correlations 
if we use the dispersion (|39| ) instead of (|^). Let us therefore use ( pOf ) to 



compute -F(n, 0) and then g{x^ 0) following Eqs. (|l^) and (pO|). The result is 
plotted and compared to the exact results for A = 2 in Fig. 8. On contrasting 
this with Fig. 1 which is based on the dispersion (|^), we see an improved 
agreement for large values of a;, and, somewhat unexpectedly, for very small 
X also. This analysis again shows the close connection between the phonon 
dispersion and the correlations in the lattice formalism. 

C. Large-x Asymptotics for g[x) 

We will now analytically evaluate g[x) in (22) for large values of x, and 
show that the order 1 and 1/A terms agree with the exact results given in 



Ref. |2^. We use the Poisson resummation formula 

OO OO 

E = dy e^'-y^ cPiy) . (40) 

n=— OO p=— OO °° 

Given a function <f){n) defined for integer values, there are clearly many ways 
to interpolate it to a function of a real variable (piy). Eq. (|iO| ) says that 
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different interpolations on the right hand side (assuming that the integrals 
and sum converge) must give the same answer for the left hand side. We will 
use this freedom to interpolate -F(n, 0) to -F(?/, 0) taking the function in (p5|), 
rather than the one in (0). This is a convenient choice because the function 
in (l25|) is smooth, unlike ( p4|) which contains oscillations. Thus, we use (^) 
with 



A 

4vr/(y) 



1/2 



exp 



4/(1/) 



(41) 



where fiy) and its first two derivatives have the form 

1 , . 1 



f{y) 



27r2 
1 



log {Cy) + 
1 



27r2y 47r4?/3 



+ 



0(: 



27r2?/2 



+ 0{ 



1 



(42) 



We now Taylor expand (j){y) around y = xpo, and do the integrals over 
y on the right hand side of (|40|). Upto order 1/A, this yields the general 
expression 

9{x) 



Pi 



f"(x) ^ / AnYfix] 



A 



A 



fix). 
;i + ^^-^) cos(27rpx) 

A 

4:1Tpf'{x) 



A 



sin(27rpx) 

(43) 



We now use the explicit forms for f{y) and its derivatives in Eq. (H2|) to 
obtain, upto order l/a;2. 



9{X) , i , r, ^^.,,l/47r2x2N-2pVA 



Pi 



2A7r2a;2 



+2 j: [c 



xe 



p=i 



1 . 
2A7r2a;2- 



cos(27rj9x) 
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2p 

sin(27rpx) 



Xttx 



(44) 

(On the right hand sides of Eqs. (^^-^), we have written x instead of xpo 
for convenience). Eq. (|4^) agrees with the asymptotic expression given in 



Ref. ||2^ upto 0{j, -^). We should make two comments here. First, the sum 



over p runs upto infinity in (|4^), while it only goes upto p = A in the exact 
result for integer A; however, the difference between the two is exponentially 
small for large A. Secondly, by going upto order 1/A^, we find terms involving 
logx which are not present in the exact result. This shows that the harmonic 
lattice approximation fails at order 1/A^, and that anharmonic terms must 
be included to recover the correct asymptotics at that order. 

The oscillatory terms in (0) imply that the Fourier transform S{q) di- 
verges at momenta q = 27rppo for all nonzero integers p satisfying < A/2. 
The divergence has a power-law form \q — 2TTppo\^^^'^^^^^^ if < A/2, and is 
logarithmic if = A/2. We indeed see a logarithmic divergence in Fig. 3 for 
A = 2, and a square-root divergence in Fig. 4 for A = 4, both at g = 27rpo- 

The asymptotic analysis can be extended to the dynamical correlations 
g{x,t). Let us define 

s = po I - vH^ 1^/^ , (45) 



and assume that s > > 1 . Then we can show that the leading term in (p!9|) is 
given by 

Fin.t) = ^ log (Cs) . (46) 
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Following the above arguments, we find that 

oo 

gix,t) = pI[i + 2J2 iCsy'^"/' cos(27rxppo) ] (47) 

p=i 

to leading order in the 1/A expansion. 

The various power-law exponents derived above are in accordance with 
arguments based on the idea of Luttinger liquids and conformal field 
theory 0. In general, if the phonons in a one- dimensional system have a 
low-momentum dispersion of the form Ug ~ vq, the arguments presented here 
show that the prefactors of the oscillatory terms in (^31) must decay with the 
exponents 27ThpQp'^/{mv). In the next section, we will study the conditions 
under which a class of long-range interactions can lead to Luttinger liquid-like 
behavior. 

III. OTHER LONG-RANGE MODELS 

We consider a one-dimensional system in which the two-body interactions 
decay as the power a of the distance. Using the lattice approximation, we find 
that the long-distance correlations are quite different depending on whether 
a > 1, q; = 1 (the Coulomb problem), or < a < 1 This directly 

follows from the phonon dispersion which has the form 



t '-^^ ■ (48) 



n=l 



(We set the density po = 1 in this section). For low momenta, g << vr, we 
see that 

Uq ~ |g| if a > 1 , 
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(a+i)/2 -f < a < 1 . (49) 



The asymptotics of the correlations g{x) are governed by the long-distance 
behavior of 

dq 



Fin) 



r - cos(gn) 1 . (50) 



This diverges as logn if a > 1 and as ^J\ogn if a = 1, but it does not 
diverge if a < 1. Correspondingly, the prefactors of the oscillatory terms in 
g{x) decay as powers of x if a > 1 (Luttinger liquid) and as exponentials of 
— ^J\ogx if a = 1, but they do not decay if a < 1 (Wigner crystal) p9|. It is 



amusing that the Coulomb system, which is physically the most interesting 
of the three cases, lies exactly in between the Luttinger liquid and Wigner 
crystal scenarios. 

IV. DISCUSSION 

We have seen that the harmonic lattice approximation is an useful tech- 
nique from which many properties of the CSM can be derived readily without 
having to solve the Schrodinger equation (|l]) in great detail. We should how- 
ever point out two problems which are either not clear or are beyond the 
reach of this approach. 

a) We have been unable to derive the one-particle Green's function, either 
bosonic or fermionic, from the phonon wave functions. This is because the 
sysmmtery (or antisymmetry) of the particle wave function does not enter 
the lattice formulation in a natural way; the lattice site and the particle 
labels are distinct from each other. This may be purely a technical difficulty, 
(b) It would be useful to systematically go to higher orders in the 1/A expan- 
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sion by starting from the original phonon dispersion (|^) and including the 
anharmonic terms perturbatively. This may provide a justification for re- 
placing the dispersion (D by either or (|39|) , and may enable us to extend 
our results to smaller values of A. 

DS thanks F. D. M. Haldane for a stimulating discussion, and the De- 
partment of Physics and Astronomy, McMaster University for its hospitality 
during the course of this work. This research was supported by the Natural 
Sciences and Engineering Research Council of Canada. 

Note Added: 

After writing this paper, we learnt that much of this work has been pub- 
lished earlier p2|, [33[]; we thank P. J. Forrester for pointing this out. Ref. |^ 



has gone much further by presenting expressions for g{x, t) at finite tempera- 
ture, for the singularities in S{q, u), and for the equal-time Green's function. 
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Figure Captions 



1. Equal-time correlation g{x)/p^ for A = 2. The solid and dotted lines 
denote the exact result and the lattice approximation respectively. 

2. Correlation g{x) / for A = 4 obtained by the lattice approximation. 

3. Static form factor S{q)/po for A = 2. The solid and dotted lines denote 
the exact result and the lattice approximation respectively. 

4. S{q)/po for A = 4 obtained by the lattice approximation. 

6. Dynamical correlation g{x,t)/pQ for A = 2 and vtpo = 0.2 obtained by 
the lattice approximation. The solid and dotted lines denote the real and 
imaginary parts respectively. 

6. Correlation g{x,t)/pQ for A = 2 and vtpQ = 2 obtained by the lattice 
approximation. The solid and dotted lines denote the real and imaginary 
parts respectively. 

7. Second moment of dynamical structure factor l2{q)m'^ / (h'^ PqQ^) for A = 2 
obtained by the lattice approximation. 

8. Equal-time correlation g{x)/pl for A = 2. The solid and dotted lines 
denote the exact result and the lattice approximation using the modified 
dispersion relation ( ^9|) respectively. 
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